Dynamics below the depinning threshold 
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We study the steady-state low-temperature dynamics of an elastic line in a disordered medium 
below the depinning threshold. Analogously to the equilibrium dynamics, in the limit T — > 0, the 
steady state is dominated by a single configuration which is occupied with probability one. We 
develop an exact algorithm to target this dominant configuration and to analyze its geometrical 
properties as a function of the driving force. The roughness exponent of the line at large scales 
is identical to the one at depinning. No length scale diverges in the steady state regime as the 
depinning threshold is approached from below. We do find, a divergent length, but it is associated 
only with the transient relaxation between metastable states. 
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The physics of elastic systems in disordered media has 
been the focus of intense theoretical and experimental 
studies. A continuing challenge has been to understand 
their response to an applied external force /. Such a 
situation is relevant for a large number of experimental 
driven systems ranging from periodic ones, such as vortex 
lattices p| and charge density waves to interfaces, such 
as domain walls in magnetic 0, Q or ferroelectric 0, 
materials, contact lines of liquid menisci Q and crack 
propagation 

A crucial feature of the zero temperature motion of 
these systems is the existence of a threshold force f c be- 
low which the system is pinned. For f > f r the sys- 
tem undergoes a depinning transition [lfj, [ll], u2, LLii LLj 
and moves with a non-zero average velocity. Fisher hrst 
viewed the depinning transition as a critical phenomenon 
A key idea is that the slow motion close to f c pro- 
ceeds by avalanches of size £, that diverge as f c is ap- 
proached. Indeed, theoretical and numerical works at 
zero temperature have shown that for / — > /+, the cor- 
relation length diverges as £ ~ (/ — / c )~ lydcp - The ex- 
ponent Vdcp is related to the roughness exponent £d cp 
of the system pinned at f c , via the scaling relation 
fdcp = 1/(2 — Cdop)- The length scale £ identifies both 
the typical size of avalanches and the crossover in the 
roughness behavior ^ij. F° r length scales below £, the 
roughness is described by the exponent Cdepj while for 
scales larger than £, the velocity introduces a time de- 
pendent noise. In this regime, for interfaces, quenched 
disorder becomes equivalent to an effective temperature, 
and the roughness exponent is equal to the thermal expo- 
nent. The interpretation of depinning as a critical phe- 
nomenon suggests that a similar diverging length scale 
would also exists below f c . Indeed, in standard critical 
phenomena, for length scales smaller than the correlation 
length the system is critical, and crosses over beyond this 
correlation length to the broken symmetry phase on one 
side of the transition and to the symmetric phase on the 
other side. An important question is thus whether an 
equivalent correlation length is observable in the limit 



Answering such a question is not an easy task since be- 
low f c the steady-state velocity vanishes in the long-time 
limit at zero temperature. This makes methods such as 
molecular dynamics simulations ill suited. Some authors 
[lri Ibil Hif have therefore studied the relaxation of a 
given initial configuration towards a final zero-velocity 
state. This analysis allows to identify a dynamical scal- 
ing characterized by a single correlation length diverging 
with an exponent ^d C p as f c is approached. However it 
has remained unclear how such T — transients relate to 
the steady state motion in the large time limit at small 
but finite temperatures. Below f c , it is thus preferable 
to keep the temperature finite, but direct Langevin dy- 
namics simulations are powerless to clarify whether the 
limit / — > /" is also characterized by a divergent length 
scale. An analytical tool for adressing this issue is the 
Functional Renormalization Group (FRG). This method 
allows to st udy, at finite temperature, the slow creep mo- 
tion [ijj, |2fJ, |2j, taking place for / <C f c . In the creep 
regime, scaling arguments rely on the physical properties 
of the system at equilibrium (/ = 0). This phenomeno- 
logical approach suggests that, for / <C f c , the macro- 
scopic forward motion is produced by activated jumps 



with typical size £r ~ / 



The exponent v eq can be 



related to the roughness by the same scaling relation, 
i/ cq = 1/(2 — C cq ), but with C eq the equilibrium roughness 
exponent. The FRG analysis 21] shows that £r play the 
role of crossover length between two roughness regimes: 
on length scales below £t, the roughness is described by 
the exponent Ccq, while for scales larger than £t it is de- 
scribed by the exponent £d cp - This regime crosses over, 
for interfaces, to a pure thermal-like behavior at very 
large length scales due to the finiteness of the velocity 
for any non-zero temperature. This FRG finding is thus 
at variance with the interpretation of the depinning tran- 
sition as a standard critical phenomenon, since it shows 
that the roughness Cdop appears at large length scales, 
and not short ones. 

In order to address the question of the behavior be- 
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FIG. 1: Dynamical phase diagram. The steady state prop- 
erties of the elastic string at T — > are determined by L opt 
(filled symbols) and £. They separate regions characterized 
by the equilibrium exponent, the depinning exponent (gray 
region), and the thermal one. The divergent length L ro i ax 
(open symbols) is associated only with transient dynamics. 
Lines are guides to the eye. 



low threshold, it is interesting to consider the small- 
temperature limit without however going all the way to 
T = 0, so that steady-state motion still exists. This 
limit is analogous to the T — > limit for the dynamics of 
a system at thermal equilibrium, where the Boltzmann 
weights impose a ground state occupation with proba- 
bility one. Occupation probabilities also exist for the 
steady-state dynamics in a finite system, but they are 
more difficult to compute. The T — > limit of this dy- 
namical regime is thus described by a single dominant 
configuration, occupied with probability one. This is 
particularly transparent for a particle hopping on a one- 
dimensional ring [231I24L |25| . 

In this paper, we construct this dominant configura- 
tions using a novel exact algorithm. We analyze its geo- 
metrical properties for driving forces between the equilib- 
rium and the depinning threshold. We find, in agreement 
with the FRG scenario for / -C f c , that the depinning 
roughness exponent describes the geometry of the line at 
large length scales for all forces. This result is inconsis- 
tent with the existence of a divergent correlation length 
for / — > /~ in the steady state regime. A divergent length 
scale exists only in the transient dynamics of relaxation 
between two metastable states. The conclusions of our 
analysis are summarized in the dynamical phase diagram 

of Fig. in 

Let us consider an elastic line on a two-dimensional 
discrete Lx M lattice with periodic boundary conditions 
in both directions (see Fig. |2J). The force / drives the 
line around the system in the direction of M. We de- 
fine elementary moves as in Ref. [2^, |2^, allowing for 
simultaneous motion of several sites. This avoids prob- 
lems of single-site dynamics proper to the elastic string 
model. Using this dynamics, any configuration relaxes to 



the nearest metastable configuration, a local minimum of 
the energy 

E = E - h ^ - f h * + v ^ w 

i 

The first term in this equation takes into account the 
short-range elastic energy and V(i, hi) is the quenched 
disorder, which we take to be Gaussian and uncorre- 
cted. The variables h gives the height of the line, as 
a function of i (see Fig. |5J. The equilibrium ground 
state (at / = 0) is easily found with a transfer matrix 
approach [3, |23] ■ The sample-dependent critical force 
f c and the zero-temperature critical configuration at f c 
are also readily determined j2(| ^ exploiting the an- 
alytic structure of the zero-temperature dynamics of a 
d-dimensional elastic interface embedded in a d + 1 ran- 
dom medium. Two properties are particularly useful for 
numerical algorithms: the 'no passing' property assures 
that, during its motion, an elastic interface can never 
miss any pinned configuration; the 'no return' property, 
which holds after a finite initial time, states that only 
forward motion takes place 30]. At non-zero tempera- 
tures, these rules obviously cease to hold. The motion is 
dominated by the activation time spent to overcome the 
energy barriers. The evaluation of these barriers is an 
NP-complete problem [3l| and the algorithms employed 
at the equilibrium are exponential in the size of the sys- 
tem. In the following we show that, in the T —> limit, 
properties analogous to the no-passing rule govern the 
time evolution of such systems. These properties allow 
to capture the steady state dynamics below the depinning 
threshold within a reasonable computation time. 

For an elastic line pinned in a metastable configura- 
tion a, the path taken at low temperature from a to 
another metastable configuration with lower energy 7 is 
the one with the lowest activation energy E^ c . The time 
spent in a thus corresponds to the Arrhcnius activation 
energy i?° sc . Once at 7, the probability to return to a 
are neglegible if T — > 0. This procedure defines a de- 
terministic coarse-grained dynamics between metastable 
states. Two theorems can be proved for this dynam- 
ics, whenever the standard no-passing rule is valid: (i) 
If there is no configuration which lowers the energy of 
a in the backward direction, the coarse-grained dynam- 
ics starting from a is always forward-directed; (ii) Let a 
be any metastable configuration escaping into a config- 
uration 7 with h 1 > h a and 7' any configuration such 
that h 1 > h a and having an energy barrier Ey c > E^ c : 

all 7' then satisfy K 1 > h 1 js^]. Direct consequences 
of these theorems are that the dynamics is periodic af- 
ter a single pass of the line around the system, and that 
the dominant configuration, the metastable configuration 
with the largest barrier, is always visited, independently 
of the initial condition. This dominant metastable con- 
figuration is the only statistically relevant configuration 
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FIG. 2: Low temperature dynamics of the driven elastic string 
below the depinning threshold: the optimal path to escape 
from a given metastable configuration a pass through a saddle 
configuration j3 that can relax deterministically to the next 
metastable configuration 7 with E n < E a . Configurations a 
and j3 differ on a length L opt ; a and 7 differ on a length L re iax- 

of the T — > steady state dynamics: Under the condi- 
tions of Arrhenius activation, the system spends much 
more time in it than in any other configuration. To con- 
struct the minimal path escaping from the configuration 
a (see Fig. 0), we enumerate all configurations dynami- 
cally connected to it, increasing gradually the maximum 
energy of configurations. This continues until a saddle 
configuration /3 is found which can relax to a metastable 
configuration 7, with E~ i < E a . Configuration a and 
/3 differ on a length L opt ; a and 7 differ on a length 
^rciax- A new path construction starts at configuration 
7. Several refinements are implemented |3^ |. During the 
construction, a very large number of configurations are 
considered, but only those with energy barriers below 
E^ c need to be evolved, fn our case, the computational 
cost increases very rapidly with L pt , but small values of 
L op t can be handled even for large system sizes L. This 
is realized in the physically interesting case close to the 
depinning threshold and the computation becomes diffi- 
cult only for small driving forces, where the saddle point 
configuration may be very different from a. We have con- 
sidered systems up to L = 128 for forces / > 0.7/ c and 
up to L = 32 for / - 0.2/ c . 

In the following, we concentrate on the disorder- 
averaged geometric properties of the dominant configura- 
tion, expressed through the structure factor S(q), which 
gives access to the roughness exponent £: 

5(g) = l^ expH^P ~ q -(^0 (2) 
3 

We also analyze the dynamical lengths L opt and L rc i a x 
defined above. In Fig. Efa), we show 5(g) for different 
forces ranging from f — (equilibrium) to / — f c (de- 
pinning). The slope of 5(g) in these two limiting cases 
corresponds to the critical exponents ( cq = 2/3 |33| and 
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FIG. 3: Steady state structure factor of the line in the T — > 
limit, averaged on 1000 samples, (a) S(q) for L = 32, M = 92 
and different forces (curves are shifted for clarity), (b) S(q) 
at f/fc = 0.8 for L = 16, 32, 64, 128 and M = L Cd <=p . 



Cdcp = 1.26 ±0.01 13], respectively. In between, we find 
two regimes. Remarkably, on large length scales (small 
g), 5(g) corresponds to Cdepj even at forces well below f c . 
In Fig. Ofb) we show that this behavior is size indepen- 
dent up to L = 128, at / = 0.8/ c . Only at small scales 
(large q) the roughness exponent is Ccq- This allows to 
define a crossover length which decreases with increasing 
/, as we can observe in Fig. [Ja). The results of Fig. |3| 
thus show that there is no divergent length scale below 
the depinning threshold in the steady state regime, unlike 
the divergent length scale that exists above the depinning 
threshold. 

In order to understand the existence of a single charac- 
teristic length which does not diverge at f c but separates 
the £eq and £dcp steady state regimes, we analyze the 
behavior of the dynamical lengths L opt and L re iax as a 
function of /. The dynamical lengths shown in Fig. ^ are 
obtained starting from the dominant configuration and 
averaged on 1000 samples. We see that L opt decreases 
with increasing / and saturates to the minimal avalanche 
size at f c . Interestingly, L rc i ax increases as (f c — f)-"^ 
close to f c and L r dax ~ L opt at low /. This behavior is 
controlled by the density of metastable states, which is 
very high at low forces, and very low near f c . The re- 
sults of Fig.^Jlead us thus to identify L opt with the steady 
state crossover length. The observed divergence in L r dax 
has its origin purely in transient dynamics and has no 
counterpart in the steady state properties. We note that 
Lrpi ax c an be related with the divergent length found in 
[id Il7l Il8l | . Our simulations allow to sketch the dynam- 
ical phase diagram of Fig. ^ This diagram is expected 
to hold for G?-dimensional manifolds in d + 1 dimensions. 
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The steady state geometrical properties of the string in 
the T — * limit are determined by the two solid lines in 
Fig. n i ptj associated with the optimal activated jump 
needed to escape from the dominant configuration, and £, 
the typical size of the depinning avalanche above thresh- 
old. L opt (f) and £(/) separate three roughness regimes: 
Coq = 2/3, Cdcp = 1-26 and Cthcrm = 0.5. The divergent 
length i rc i ax (dashed line) does not affect steady state 
properties: it describes transient processes that depend 
on the distance between successive metastable states, but 
it does not describe the properties of the dominant ones. 
As indicated in Fig. it is plausible to connect contin- 
uously L opt with the thermal nucleus size £t defined for 
/ <C f c - However, we are not able to verify with our data 
the predicted power law divergence £t ~ /~ !y ° q - The sce- 
nario thus proposed is consistent with the FRG predic- 
tion |2lJ of the existence of Cdep at scales larger than £t 
and inconsistent with the existence of a critical region ex- 
tending both below and above / c , as the standard critical 
phenomena interpretation would have suggested. 

Although our computational approach strictly holds 
only in the T — > limit, it yields the scenario of Fig. ^ 
which offers a solid framework for discussing finite tem- 
perature effects. At finite temperature the depinning 
transition is rounded. Since the velocity is finite for all 
forces, the regime with the thermal roughness Ctherm ex- 
ists at the largest length scales |21|. For / < f c we ex- 
pect two crossovers separating the regimes Ccq, Cdcp and 
Cthcrm, as we increase the length scale. For temperature 
comparable to the strength of the disorder at the very 
small length scales an additional thermal regime appears 
and crosses over to Ccq or directly to Cdcp- Such a sce- 
nario is in part supported by recent Langevin dynamics 
simulations (2^, which show, at small forces, a steady 
state motion characterized by a roughness exponent big- 
ger than Ceq at larges length scales and a roughness Ctherm 
at the smallest length scales. 

Finally, we expect that an experimental verification of 
our results is possible using, for instance, imaging tech- 
niques for magnetic 0,0] or electric 0,E| domain walls in 
thin films: L opt could be extracted from the analysis of a 
spatial correlation function or by measuring the domain 
wall speed at low temperatures , si nce Lopt also controls 
the thermally activated motion |32j ; £ re iax could be mea- 
sured by comparing consecutive (long-lived) metastable 
states when / is close to f c , since then L rc iax controls the 
distances between successive metastable states. 
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